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Abstract 

The Coastal Atmosphere and Sea Time Series (CoASTS) Project, aimed at supporting ocean color research 
and applications, from 1995 up to the time of publication of this document, has ensured the collection of a 
comprehensive atmospheric and marine data set from an oceanographic tower located in the northern Adriatic 
Sea. The instruments and the measurement methodologies used to gather quantities relevant for bio-optical 
modeling and for the calibration and validation of ocean color sensors, are described. Particular emphasis 
is placed on four items: 1) the evaluation of perturbation effects in radiometric data (i.e., tower-shading, 

instrument self-shading, and bottom effects); 2) the intercomparison of seawater absorption coefficients from 
in situ measurements and from laboratory spectrometric analysis on discrete samples; 3) the intercomparison 
of two filter techniques for in vivo measurement of particulate absorption coefficients; and 4) the analysis 
of repeatability and reproducibility of the most relevant laboratory measurements carried out on seawater 
samples (i.e., particulate and yellow substance absorption coefficients, and pigment and total suspended matter 
concentrations). Sample data are also presented and discussed to illustrate the typical features characterizing 
the CoASTS measurement site in view of supporting the suitability of the CoASTS data set for bio-optical 
modeling and ocean color calibration and validation. 


1. INTRODUCTION 

Optical measurements from space enable estimates of 
the concentration of materials suspended or dissolved in 
seawater (i.e., pigment, sediment, and colored dissolved 
organic matter) which are of great relevance in environ- 
mental and climate-related studies. Because of this, a 
number of new advanced ocean color sensors were designed 
to support oceanographic studies and applications includ- 
ing the Modular Opto-electronic Sensor (MOS), the Ocean 
Color and Temperature Scanner (OCTS), the Polarization 
and Directionality of the Earth s Reflectance (POLDER) 
sensor, the Sea-viewing Wide Field-of-view Sensor (Sea- 
WiFS), the Moderate Resolution Imaging Spectroradiome- 
ter (MODIS), and the Medium Resolution Imaging Spec- 
trometer (MERIS). All of these sensors were successfully 
launched and have contributed significantly to the general 
problem of inverting optical measurements to derive con- 
centration estimates of biogeochemical parameters; several 
continue to provide regular coverage of the global bio- 
sphere. 

To ensure a better exploitation of the data supplied by 
spaceborne sensors, national and international calibration 
and validation projects have been started. Their major ob- 
jectives are the development of the bio-optical algorithms 
required for extracting quantitative information from space 
data; the validation of products obtained from satellite im- 
agery; and the indirect absolute calibration (i.e., vicarious 
ciibration) of the radiometers in space. All of these activ- 
ities require comprehensive in situ atmospheric and ma- 
rine measurements. In agreement with this requirement, 
the Coastal Atmosphere and Sea Time Series (CoASTS) 
Project was set up at the Joint Research Centre (JRC) 
in Ispra, Italy to support bio-optical modeling and ocean 
color calibration and validation exercises at a coastal site 
in the northern Adriatic Sea. 


In order to provide a comprehensive overview of the 
CoASTS measurement activities, this report has three pri- 
mary objectives: 

1) Present the instrumentation and the measurement 
methodologies, 

2) Assess the accuracy of the most relevant measure- 
ments needed for bio-optical modeling as well as 
ocean color vicarious calibration and algorithm val- 
idation activities; and 

3) Discuss sample data that display the relevant fea- 
tures of the measurement site. 

2. BACKGROUND 

Ocean color calibration and validation activities carried 
out in the late 1970s and early 1980s for the Coastal Zone 
Color Scanner (CZCS), relied exclusively on a small num- 
ber of ship-based campaigns for ground truth observations. 
More recent calibration and validation activities for MOS, 
OCTS, SeaWiFS, and MODIS make use of moored buoys, 
in addition to an extensive suite of oceanographic research 
cruises. Examples of radiometric buoy systems that were 
successfully deployed are as follows: 

■ The Yamato Bank Optical Mooring (YBOM), de- 
veloped by the National Space Development Agency 
of Japan (NASDA) for OCTS calibration and vali- 
dation (Kishino et al. 1997); 

■ The Marine Optical Buoy (MOBY), developed by 
the National Aeronautics and Space Administration 
(NASA) and the National Oceanographic and At- 
mospheric Administration (NOAA) to support Sea- 
WiFS and MODIS calibration and validation activ- 
ities (Clark et al. 1997); and 
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■ The Plymouth Marine Bio-Optical Data Buoy 
(PlyMBODy), developed by the Plymouth Marine 
Laboratory (PML) in the United Kingdom for Sea- 
WiFS ocean color calibration and validation (Pink- 
erton and Aiken 1999). 

More recently, a new European radiometric buoy is being 
deployed as part of the Bouee pour V acquisition de Series 
Optiques a Long Terme] (BOUSSOLE) Program (Antoine 
and Guevel 2000). 

YBOM was equipped with a fluorometer, two pairs of 
upwelling radiance and downward irradiance sensors at 1.2 
and 6.5 m depths providing spectral data from 400-800 nm, 
and an above-water irradiance reference sensor providing 
data at 860 nm. YBOM was first deployed and tested 
in the latter part of 1993, redeployed almost two years 
later for final testing, and operationally used during the 
Advanced Earth Observing Satellite (ADEOS-1) mission 
(Shimada et al. 1998). 

The MOBY program became operational in July 1997 
with three buoy systems supported by a deployment and 
refurbishment facility established at the University of Ha- 
waii Marine Center in Honolulu. The buoy systems are de- 
ployed on a rotational basis off the coast of Lanai, Hawaii. 
The site is in the open ocean (low chlorophyll a concentra- 
tion) with a constant aerosol type, and ensures availabil- 
ity of high resolution data throughout the SeaWiFS mis- 
sion. MOBY is equipped with an above- water irradiance 
reference sensor and three pairs of downward irradiance 
and upwelling radiance collectors located at 2, 5, and 9 m 
depths, all providing data from 380-900 nm using custom- 
built hyperspectral sensors. The light sensors are mounted 
on arms to minimize the self-shading effect of the buoy, and 
orientation sensors allow the exclusion of data when the 
optical sensors are behind the buoy or excessively tilted. 
Mitigation of biofouling is accomplished with an aggres- 
sive campaign of diver visits to keep the optical windows 
clean. Coincident casts with an optical profiler are used to 
monitor data degradation between cleaning sessions and to 
independently demonstrate the baseline capabilities of the 
buoy data when the optical surfaces are not contaminated. 

PlyMBODy was equipped with a conductivity, temper- 
ature, and depth (CTD) sensor; a fluorometer; an above- 
water irradiance reference sensor; and a pair of upwelling 
radiance sensors at 0.7 and 2.6 m depths. The radiometers 
were manufactured by Satlantic, Inc. (Halifax, Canada), 
and provided data in seven spectral bands from 412- 
780 nm. PlyMBODy was deployed off Plymouth in the 
English Channel, and was designed for a three-year life- 
time from 1997-1999 (assuming a 9-month deployment 
from April-December). The full-scale beta version was 
deployed on site during a 3-month period (20 July to 20 
October 2000), which allowed it to encounter a variety of 
meteorological situations. The buoy was equipped with 


t Buoy for the acquisition of a long-term optical series. 


two inclinometers, a pressure sensor, an ARGOS beacon, 
and a flashing light. The goal of this deployment was to 
establish the capabilities of the buoy as well as to identify 
possible problems necessitating modifications. 

The BOUSSOLE buoy will be deployed in the Ligurian 
Sea (between France and Corsica) in early 2002 as part 
of the MERIS calibration and validation activities being 
developed at Laboratoire d’Oceanographie de Villefranche 
(LOVf), located in Villefranche-sur-Mer, France. In terms 
of the optical measurements, the system consists of Sat- 
lantic radiometers measuring the solar irradiance (at 4.5 m 
above the surface), plus the downward irradiance, upward 
irradiance, and upwelling radiance at 4 and 9 m depths. 
The two-axis tilt, atmospheric pressure, and buoy orien- 
tation with respect to the sun, are also recorded. In ad- 
dition, two fluorometers are mounted at 4 and 9 m, plus 
a transmissometer and a CTD are deployed at 9 m. Two 
sets of in-water optical sensors are used to minimize bio- 
fouling. While one set is in use, the other is characterized 
while fouled, cleaned, and then characterized again before 
redeployment — approximately one month after extraction 
by divers. 

The use of buoys or ships has advantages and draw- 
backs. Moored buoys can easily ensure an almost continu- 
ous collection of a restricted set of parameters at a specific 
site. Alternatively, ship-based measurements can ensure a 
wide spatial data collection with a comprehensive optical 
and biogeochemical characterization of each sampling site, 
but with poor temporal resolution. 

Recently, the use of oceanographic towers, as an alter- 
native to buoy or shipboard measurements, was investi- 
gated (Zibordi et al. 1995 and Kearns et al. 1996). Com- 
pared to buoys, the use of offshore towers as mooring sites 
for continuous measurements has the disadvantage of in- 
creased contamination of the optical data with superstruc- 
ture perturbations. The use of oceanographic towers as 
logistic platforms for comprehensive optical and biogeo- 
chemical measurements, when compared to ships, has the 
disadvantage of not permitting any spatial observation ca- 
pability. 

Despite the limitations, oceanographic towers offer the 
unique opportunity of a very stable measurement platform 
enabling easy and complete control of the deployment ge- 
ometry (Hooker et al. 1999 and 2002). The latter feature 
may ensure the deployment of optical instruments with 
virtually no tilt and the exact identification of the solar il- 
lumination geometry needed for an accurate removal of su- 
perstructure shading effects (Zibordi et al. 1999). Taking 
advantage of these deployment conditions, the CoASTS 
Project has been relying on the use of the Acqua Alta 
Oceanographic Tower (AAOT) in the northern Adriatic 
Sea for the collection of a comprehensive atmospheric and 


i Formerly the Laboratoire de Physique et Chimie Marines 
(LPCM). 
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Fig. 1 . The number of optical stations, N, performed in each CoASTS campaign (numbered sequentially 
above the bars) for the period October 1995'-December 1998. 


marine coastal data set. Periodic short deployments, us- 
ing a modular in-water system that permits the removal of 
the light sensors in between measurement campaigns, will 
completely eliminate any biofouling of the in- water optical 
sensors. 

3. COASTS OBJECTIVES 

The CoASTS Project, started in 1995 and still contin- 
uing, has three aims: 

1) The collection of the time series of atmospheric and 
marine data in the northern Adriatic Sea from a 
single measurement site for ocean color research and 
applications; 

2) The analysis of the acquired data for the develop- 
ment of site- and seasonal-specific marine and at- 
mospheric algorithms for ocean color data exploita- 
tion; and 

3) The use of the acquired data in calibration and val- 
idation activities. 

Recent studies (Berthon et al. 2000, Hooker et al. 2000a, 
Zibordi and Berthon 2001, Berthon et al. 2002, Sturm and 
Zibordi 2002, and Zibordi et al. 2002) have shown the capa- 
bility of the project in reaching the former objectives and 
has clearly established the benefits offered by an ongoing 
time series of measurements (i.e., the continuous increase 
of the statistical significance of the data set for bio-optical 
modeling, and the expansion of the number of matchups 
for validation or vicarious calibration exercises). 

4. MEASUREMENT PLAN 

The CoASTS field campaigns, performed onboard the 
AAOT located 8 miles southeast off the Venice Lagoon 
(12.510°E,45.310°N), take place for a few days every 2 
4 weeks. The frequency of campaigns is generally higher 
when biological phenomena are more pronounced (i.e., in 
the boreal spring and fall), and lower when environmental 


effects decrease the likelihood of good optical data or safe 
embarkation onto the tower (i.e., during the boreal winter). 
Each comprehensive data collection — called a “station” 
includes the in-water optical and hydrographic profiles, 
seawater samples at different depths (surface, 8 m, and 
14m), atmospheric optical measurements, meteorological 
data, and visual observations of cloud cover and sea state. 

For each station, data and samples are generally col- 
lected in less than 30 min. Treatment of seawater samples 
(i.e., filtration and storage of samples in liquid nitrogen or 
in a refrigerator) generally requires 1 2 extra hours. Labo- 
ratory sample analysis is then performed within a few days 
of concluding the measurement campaign. 

Figure 1 shows the number of campaigns and the num- 
ber of stations per campaign starting in the latter part of 
1995 and extending for three years. The number of sta- 
tions per campaign, varying from a minimum of 3 up to a 
maximum of 10, is mostly determined by the meteorologi- 
cal conditions that may restrict the access to the tower or 
largely reduce the number of possible stations performed 
on ideal illumination conditions (i.e., when the sun is not 
covered by clouds or with an overcast sky ) . 

5. SITE CHARACTERISTICS 

The CoASTS site is located in a frontal region that can 
be characterized by Case-1 or Case-2 water types (Berthon 
et al. 2002). The Case-2 water characteristics are mostly 
determined by the input from the northern rivers (i.e., Pi- 
ave, Livenza, and Tagliamento) . The aerosol type, occa- 
sionally maritime, is mostly continental and is determined 
by atmospheric input from the nearby Po River valley. 
These characteristics of the site, represent most of the 
northern Adriatic Sea region. 

The horizontal homogeneity in the vicinity of the mea- 
surement site, evaluated through SeaWiFS images, has dis- 
played extremely different situations. The standard devi- 
ation of the satellite- derived normalized water-leaving ra- 
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Table 1 . The typical average values of quantities characterizing the measurement site. The ± values indicate 
standard deviations. 


Quantity 

C a [//gL *] 

Ctsm [mgL ! ] 

KE d (490) [m- 1 ] 

o vs (400) [m 

t^(555) 

Average 

1.3 ±1.1 

1.1 ± 0.7 

0.21 ±0.09 

0.15 ±0.06 

0.14 ± 0.06 


diance at 555 nm, Lwn{ 555), computed from a 5x5 ma- 
trix of image pixels centered at the AAOT site, has shown 
values ranging from 0.02-0.66 (Sturm and Zibordi 2002) 
in the range of 0.5-2. 5mWcm~ 2 /im -1 sr -1 . This result 
indicates that at the measurement site, there are time 
periods characterized by high spatial inhomogeneity and 
other time periods characterized by high spatial unifor- 
mity. When the latter is the case, the AAOT in situ data 
can be assumed to be representative of the area surround- 
ing the measurement site. 

Table 1 presents the typical values of the most rele- 
vant quantities characterizing the CoASTS measurement 
site: the chlorophyll a concentration, C a \ the total sus- 
pended matter concentration, Ctsm; the diffuse attenua- 
tion coefficient for downwelling irradiance at 490 nm, Ks d ; 
the colored dissolved organic matter (CDOM) or “yellow 
substance” absorption coefficient, a ys , at 400 nm; and the 
aerosol optical thickness, r a , at 555 nm. 

6. THE AAOT 

The AAOT, built in 1970, is owned by the Istituto per 
lo Studio della Dinamica delle Grandi Masse f (ISDGM) 
of the Italian National Research Council (CNR) in Venice. 
The AAOT, mounted on four grounded pillars, has four 
levels. Each level is 7.2 mx 5.2 m in size with the exception 
of the lowest level which is 5.2 m 2 . Measurement activities 
are generally carried out on the second and fourth levels. 

On the second level, at approximately 7 m above the sea 
surface, a portable laboratory provides space for treating 
water samples and for data logging. On the same level, an 
open grid platform 3.5 m wide, extends 6.5 m over the sea 
facing southeast; it also provides mounting points for the 
instruments to be deployed into the sea. 

On the fourth level, at approximately 12 m above the 
sea surface, the meteorological and atmospheric optical in- 
struments are installed. Electrical power is provided by 
diesel-powered generators and by lead-acid batteries. The 
generators provide electrical power during field campaigns 
when the AAOT is being used, and the batteries ensure 
power for the instruments used in continuous measurement 
programs. 

7. THE MEASUREMENTS 

The CoASTS data collection activities are split into 
field and laboratory work. The field measurements taken 
are composed of the following: 


f The Institute for the Study of Dynamics of Large Masses. 


a) In-water, vertical profiles of spectral (z indicates 
depth and A indicates wavelength) upwelling nadir 
radiance, L u (z, A), downward irradiance, E d (z, A), 
and upward irradiance, E u (z , A), measured with the 
Wire-Stabilized Profiling Environmental Radiome- 
ter (WiSPER) composed of Satlantic components; 

b) Above-water total solar irradiance, E d ( 0 + ,A), and 
diffuse sky irradiance, Ei(0 + ,A), recorded with a 
Yankee Environmental Systems (YES), Inc. (Turn- 
ers Falls, Massachusetts) MFR-6 rotating shadow 
band radiometer (z = 0 + indicates an above- water 
measurement); 

c) Direct solar irradiance, E( A), and diffuse sky radi- 
ance, Li(0, 4 i>, A), at different zenith, 0 , and azimuth, 
0, angles in the almucantar and sun planes, taken 
with a CIMEL Electronique (Paris, France) CE-318 
sun photometer; 

d) Profiles of seawater beam attenuation, c(z, A), and 
absorption, a(z, A), coefficients, taken with a West- 
ern Environmental Technology Laboratories ( WET- 
Labs), Inc. (Philomath, Oregon) AC-9; and 

e) Ancillary field data (i.e., profiles of seawater tem- 
perature, T w (z), and salinity, S^z); Secchi disk 
depth, SD; tide level, TJ; atmospheric pressure, P a \ 
relative humidity, RH ; air temperature, T a ; wind 
speed, W s , and direction, W d \ cloud cover, OC\ and 
sea state, M). 

The laboratory measurements, which provide data com- 
plementary to the field measurements, are: 

■ In vivo particulate absorption coefficient a p (A) sep- 
arated into its two components a p h(A) for the pig- 
mented and Qdp(A) for the nonpigmented matter, 
obtained from spectrometric analysis of particles re- 
tained on filters; 

■ CDOM absorption coefficient a ys (A) obtained from 
spectrometric analysis of filtered seawater; 

■ Pigment concentration from high performance liq- 
uid chromatography (HPLC) analysis; 

■ Total suspended matter concentration, Ctsm, de- 
termined with the dry weighting technique; and 

■ Particle size distribution (PSD) spectra from Coul- 
ter Counter analysis. 

In the collection of the CoASTS measurements, great 
effort was placed in ensuring consistency of data in a vari- 
ety of ways: 

1) Using assessed protocols for data collection, calibra- 
tion, and preprocessing; 


4 




G. Zibordi, J.-F. Berthon, J.P. Doyle, S. Grossi, D. van der Linde, C. Targa, and L. Alberotanza 


2) Choosing data collection frequencies aimed at re- 
solving seasonal variations in the observed quanti- 
ties; 

3) Assuring uniformity in the accuracy of the data over 
time; and 

4) Attempting regular collection of each quantity in- 
cluded in the time series. 

Most of the former data consistency elements, are eas- 
ily achieved in a short-term measurement program, but 
over long periods (years) they can be affected by a number 
of elements: 

• Old instruments being replaced by new instruments 
that have slightly different features (i.e., spectral 
response, sensitivity, or calibration uncertainties); 

• Instrument failures temporarily interrupting specific 
measurements; 

• Occasional harsh environmental conditions prevent- 
ing regular data collection; and 

• Revisions to the existing measurement protocols be- 
cause of advances in the state of the art. 

8. INSTRUMENTS AND METHODS 

With the goal of highlighting the consistency of the 
CoASTS Project resulting from a multiyear data collec- 
tion program, the following subsections describe the mea- 
surement methods and present the different quality control 
elements. The latter include, but are not limited to, abso- 
lute calibration accuracy, repeatability and reproducibility 
analysis, and intercomparison exercises. Particular impor- 
tance is given to uncertainty analyses of radiometric ma- 
rine quantities relevant for both bio-optical modeling and 
calibration and validation activities, i.e., L u (z , A), E d (z , A), 
and E u (z , A). 

8.1 Optical Measurements 

In the following subsections, the instruments for in- 
water radiometric data collection are described together 
with the measurements and absolute calibration schemes, 
the techniques for subsurface values computation, and re- 
moval of tower-shading, self-shading, and bottom pertur- 
bations. An attempt is also made to quantify environ- 
mental effects (i.e., sea state, changes in seawater optical 
properties, and illumination conditions) on subsurface ra- 
diometric quantities. 

8.1.1 Instruments and Measurement Methods 

Measurements of L u (z, A), E d (z, A), and E n {z,\) are 
obtained in seven spectral bands with a 10 nm bandwidth. 
The WiSPER system is installed on a metal frame which 
moves along two vertical wires fixed between the tower 
and the sea bottom (Fig. 2). The wires provide a high 
degree of stabilization of the WiSPER instrument system 
and prevent almost all movement out of the vertical plane. 
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Fig. 2. A schematic of the AAOT showing the rele- 
vant items to the CoASTS program: a) the deploy- 
ment platform for atmospheric instruments, b) the 
laboratory for water filtration and data logging, c) 
the deployment platform for in-water instruments, 
and d) the rig for in-water profiling instrumenta- 
tion. The inset panel displays the AAOT top view. 

The WiSPER radiometers are Satlantic Ocean Color 
Radiance and Irradiance series-200 sensors, OCR-200 and 
OCI-200, respectively. One OCR-200 is used to measure 
L u {z, A), and one OCI‘200 is used to measure E d (z, A) 
and E u (z, A), where E u {z, A) is measured by rotating the 
bracket holding the OCI-200 for E d (z, A) 180°. The light 
sensors are mounted on an extension boom, which places 
them 1 m away from the main part of the frame and ap- 
proximately 7.5 m from the tower legs. Power, analog-to- 
digital (A/D) conversion, and telemetry of the light sensor 
data are provided by a Satlantic (16-bit) module called a 
DATA- 100. 

The effective signal-to-noise ratio (SNR) for L u (z, A) 
in the spectral range 412-555 nm, is generally higher than 
4xl0 3 , and more than 2xl0 3 at 665 nm; the effective SNR 
for E d {z,\) i s generally greater than 10 4 . Because the 
same OCI-200 is used for both E d {z } A) and E u {z,\), the 
effective SNR for E u {z, A), only varies from 5-10xl0 2 in 
the spectral range 412-555 nm, and is about 1.5xl0 2 at 
665 and 683 nm. 

WiSPER, together with the AC-9 installed on the same 
rig to provide simultaneous measurements, is raised and 
lowered from the southeastern side of the tower by an elec- 
trical winch. The speed of the winch is approximately 
0.1ms -1 , which was chosen to satisfy the AC-9 deploy- 
ment requirements (i.e., to have a water flux in the AC-9 
measurement chambers ensuring the characterization of 
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Table 2 The nominal center wavelengths (in nanometers) for the channels of the in situ optical instruments 
used in the CoASTS Project. 


Instrument and Channel 

1 

2 

3 

4 

5 

6 

7 

8 

9 

WiSPERf L u (z, A) 

412 

443 

490 

510 

555 

665 

683 



WiSPERj E u (z, A) and Ed(z , A) 

412 

443 

490 

510 

555 

665 

683 



AC-9 a(z, A) and c(z, A) 

412 

440 

488 

510 

555 

630 

650 

676 

715 

MFR-6 E d { 0 + , A) and £ t (0 + , A) 

415 

500 

610 

665 

862 

960 




CE-318 Li(0, <f>, A) 



440 


675 

870 

946 

1020 


CE-318 E( A) 

+ A C + 1 — „ Xfimc . 

340 

380 

440 

501 

675 

870 

1020 




j A Satlantic MVDS reference to measure E d ( 0 + ,A) with the same nominal center wavelengths was added in the second 


seawater optical properties with a depth resolution better 
than 0.2 m). 

During each measurement station, three sets of up-and- 
down profiles are made: the first and the third to collect 
L u (z, A) and Ed(z, A), and the second to collect L u (z, A) 
and E u (z, A). Significant variations between quantities de- 
rived from the first and the third profile are used as an 
indicator of unstable environmental conditions. The sec- 
ond profile is primarily used to determine the Q-factor at 
nadir viewing: 


Qn{z, A) 


E u jz, A) 

Ms, A) - 


( 1 ) 


All of the collected data are automatically binned within 
the DATA- 100 with a 0.2 m vertical resolution. 

The major drawback of the applied methodology is 
the lack of continuous above-water irradiance measure- 
ments for minimizing the effects of illumination changes 
during casts. To overcome the former limitation, a Sat- 
lantic Multichannel Visible Detector System (MVDS) so- 
lar reference system was added to the instrumental set in 
the latter half of 1998. The MVDS provides continuous 
£d(0 + . A) measurements at the same time of in- water pro- 
filing. Moreover, the WiSPER system was upgraded by 
adding a second OCI-200 to permit simultaneous measure- 
ments of L u (z> A), Ed(z, A), and E u (z, A). A summary of 
the optical instruments used in the CoASTS Project is 
given in Table 2. 

The following work only deals with data collected up 
to the latter half of 1998 with the two-sensor WiSPER 
configuration, but most of the conclusions are applicable 
to the data collected with the three-sensor configuration 
without any loss in generality. 


8.1.2 Absolute Calibration 

The absolute calibration of the considered radiometric 
quantities is obtained from 

3(2, A) = [V-( 2 ,A) - D t { A)] c 5 ( A)/ 5 (A), (2) 

where the “hat” notation ( A ) indicates a measured quan- 
tity. The 3 symbol (which is used both as a variable and a 


subscript) indicates the calibrated radiometric quantities! 
M*.A), Ed(z ) A), and E u (z, A) (upwelled radiance, down- 
ward irradiance, and upward irradiance, respectively) not 
yet corrected for any perturbation effects (i.e., self-shading, 
tower-shading, and bottom). The digital voltage (in counts) 
at center wavelength A and depth z is V$(z, A), D$( A) is 
the average dark value (taken at the end of each measure- 
ment station by putting caps on the radiometers), C^(A) is 
the absolute calibration coefficient, and i^(A) is the immer- 
sion factor (i.e., the coefficient accounting for the change 
in sensor responsivity when the in-air calibration is applied 
to in-water measurements). 

The absolute radiometric calibration coefficients — 
an< ^ ^E d (A) [which are identical in the two-sensors 
WiSPER configuration, because of the use of the same 
radiometers for Ed(z, A) and E u (z, A) measurements] — are 
determined using a 1,000 W FEL lamp traceable to the Na- 
tional Institute of Standards and Technology (NIST). The 
calibration coefficients (A) Eire determined making use 
of an FEL lamp and a Labsphere (Sutton, New Hampshire) 
18 in Spectralon 99% reflectance plaque. In order to track 
any temporal change in the sensitivity of a light sensor, the 
calibration coefficients are generally determined every 3-6 
measurement campaigns. When sensitivity changes larger 
than 1.5% are detected, a linear variation of the sensor sen- 
sitivity is assumed from one calibration date to the next 
as a function of the number of measurement campaigns 
falling within the two dates. 

The uncertainty in the absolute calibration coefficients 
Cj(A) is ±2.1% for radiance (mostly resulting from the 
quadrature sum of the lamp calibration and plaque re- 
flectance uncertainties, both taken equal to ±1.5%) and 
±1.5% for irradiance (mostly due to the uncertainty in 
the lamp calibration values). These estimated values are 
confirmed by the detailed analyses of radiometric absolute 
calibration uncertainties presented by Hooker et al. (2002). 

The immersion factors used in the calibration relation- 
ship are provided by the instrument manufacturer (in this 


t For the CoASTS data set, radiance measurements are in 

units of mW cm 2 /xm~ 1 sr“ 1 and irradiance measurements 
are in units of mWcm -2 /xm _1 . 


6 




G. Zibordi, J.-F. Berthon, J.P. Doyle, S. Grossi, D. van der Linde, C. Targa, and L. Alberotanza 


case, Satlantic). The immersion factors for radiance sen- 
sors were computed theoretically following Mueller and 
Austin (1995), while the irradiance immersion factors were 
determined experimentally (S. McLean, pers. comm.). Be- 
cause of the lack of independent and assessed data, the 
accuracy of immersion factors was assumed to be ±1% for 
radiance (mostly due to the uncertainty in the refractive 
index of the optical window) and ±3% for irradiance [us- 
ing uncertainties estimated for cosine collectors made of 
the same material as the OCL200 (Mueller 1995)]. 

In addition to uncertainties for C^(A) and I^( A), the 
non-ideal cosine response of the irradiance collectors could 
be a further source of uncertainty in irradiance measure- 
ments. In the specific case of the OCI-200 radiometers, 
the non-cosine response may cause slightly different un- 
certainties at different wavelengths because of the use of 
separate collectors for each radiometric channel. The Sea- 
WiFS measurement protocol (Mueller and Austin 1995) 
requires the uncertainty in the irradiance measurement of 
a collimated light to be better than 2% for angles from 
0-65° with respect to the optical axis, and within 10% for 
angles from 65—90°. Because of the unavailability of pub- 
lished data on the accuracy of the cosine response for Sat- 
lantic OCT 200 radiometers, a ±2% uncertainty is assumed 
at all wavelengths. 

The uncertainties for C^(A), I^( A), and the cosine re- 
sponse of irradiance sensors, lead to a quadrature sum of 
±2.3% and ±3.9% for the overall radiance and irradiance 
uncertainty values, respectively. 

8.1.3 Computation of Subsurface Values 

A subsurface valuef 3(0“, A) [i.e., L u (0 ,A), Ed { 0 , A) 
and £ U (0“,A)] is the exponent of the intercept from the 
least-squares linear regression of ln(3(^,A)) versus depth 
within a near-surface interval zq <z<z\, where ln(3(^ A)) 
is either ln(L u (z, A)), In (Ed(z, A)), or In (E u (z, A)). Depths 
zo and z\ are iteratively chosen (D’Alimonte et al. 2001) 
to satisfy the requirement of a linear decay in ln3(^. A) as 
a function of z (generally for the CoASTS stations 0.3 < 
Zq < 1 m and 2.5 < z x < 4.5m). The negative value of the 
slope of the regression fit is the diffuse attenuation coeffi- 
cient iCj(A) [i.e., Ke u (X) and Kg d ( A)]. 

The subsurface values are corrected for perturbation 
effects according to 

L u (0',A) = Lu(0 _ ,A)r ? F 5 (A)r7ff(A)r ? f B (A), (3) 

E d ( 0-,A) = Ed(0-,\)T%*(A), (4) 

and 

£ U (0',A) = £ u (0',A)r ? ^(A)7 7 |®(A)7/| u B (A), (5) 


f A depth immediately below the sea surface is indicated by 
z = 0‘. 


where L u ( 0 _ , A), E d { O', A), and E u ( O', A) are the corrected 
values, and r)7 s (\), 7/| s (A), and r?? E ( A) are the tower- 
shading, self-shading, and bottom effects correction fac- 
tors, respectively. 

The application of the correction factors in (3)-(5) as- 
sumes the sources of the perturbations are all indepen- 
dent. The accuracy of the corrections is also based on 
the assumptions that the tower-shading and self-shading 
effects are depth independent, and the bottom effects are 
weak. Because L u ( 0”,A), Ed ( 0 , A), an( 3 E u (0 , A) are all 
computed from subsurface data (generally within the 0.3- 
4.5 m depth interval) and bottom effects are, in most of the 
cases, less than 2% of the signal at z = 0 , the correction 
formulations given in (3)-(5) can be applied to most of the 
CoASTS stations. 

The factor rj ^ associated with each perturbation effect 
is defined as 


^(A) 


3(0", A) 
3(0", A) ’ 


(6) 


where 3(0", A) is the measured radiometric quantity and 
3(0", A) is the corrected value [i.e., L u (0",A), E d ( 0",A), 
and E u (0~ , A)]. 

The uncertainties related to perturbations discussed in 
this document, make reference to data taken during so- 
called clear-sun, clear-sky, and almost-stable environmen- 
tal conditions. During clear-sun conditions, there are no 
clouds in front of, or in near proximity of, the solar disc. 
During clear-sky conditions, the sky is cloudless (except 
perhaps at the horizon). Clear-sun and clear-sky condi- 
tions are discriminated by a threshold for the ratio be- 
tween above- water diffuse and direct irradiance at 412 nm, 
r*(412), obtained from the MFR-6 data measured at the 
same time as the in-water optical profiles; the threshold 
used for data processing purposes is n( 412) < 2.0. Almost- 
stable environmental conditions are identified by differ- 
ences lower than 10% between the L u ( 0 ,\)/E d (0 , A) ra- 
tios computed from successive radiometric profiles (i.e., the 
first and the third) taken in each measurement station. 

Figure 3 shows the temporal evolution of the total cor- 
rection factors tj 7 s (^) f?^ S (A) applied to L u ( 0 ,A) 

and the frequency distribution of the related percentage 
values 100[t£S(A)^f(A)i7*f(A) - l] at 412, 490, 555, 
and 665 nm (where the latter wavelengths were chosen to 
ensure that the correction factors are represented over the 
visible spectrum). The total correction factors exhibit sig- 
nificant spectral dependence and very large variations at all 
wavelengths, with minimum values ranging from 0.91-1.10 
at 555 nm, and maximum values ranging from 1.06-1.16 
at 665 nm. The frequency distribution of percentage cor- 
rections shows average values of 9.2, 5.5, 3.2, and 9.3% at 
412, 490, 555, and 665 nm, respectively. The former results 
highlight the spectral and seasonal dependence of pertur- 
bation effects which, if not removed, can induce significant 
inconsistencies in the time series of L u ( 0 , A), E d ( 0 ,A), 
and E u ( O', A). 
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Fig. 3. The total correction factors for tower-shading, self-shading, and bottom effects vf S ( A) vf E ( A) 

apphed to upwelhng radiance as a function of the station date at A = 412, 490, 555, 665 nm, and" the related 
percentage frequency distribution [ty|?(A) r;| 5 (A) t?| £ (A) - l] 100. The plus symbol (+) indicates single 
station values, while the dot symbol (•) indicates mean campaign values. The average is indicated as av, and 
c 7 is the standard deviation. The number of points is n — 215. 
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Table 3. The reference values used for computing the tower-shading correction factors at each OC 1-200 and 
OCR-200 nominal center wavelength (given in units of nanometers) for L u (0 ,A), E u ( 0 ,*), and ^d{ l J , )■ 


Index 

412 

443 

490 

510 

555 

665 

a [m 

x ] 

0.02 

0.05 

0.10 

0.30 

0.50 



0.50 

0.70 

0.80 

0.85 

0.90 

0.95 

6 0 [°] 


25 

30 

40 

50 

60 

70 

00 [°] 


90 

135 

180 

225 

270 


r a 


0.00 

0.05 

0.10 

0.50 

1.00 

— + oc 


Table 4. The OCI-200 and OCR-200 wavelength-specific values used for computing tower-shading correction 
factors (the wavelengths are given in units of nanometers). 

' — 1 “ I 


iatlUl O y tilV. LA. » N. 

Parameter 

O* 

412 

443 

490 

510 

555 

665 

PB 

To 

TR 

1.3xl0 -2 1 

6.6x 10 -4 
0.32 

1.8x 10 -2 
l.OxlO -3 
0.24 

3.1xl0 -2 

7.4xl0 -3 

0.16 

4.0 x 10 -2 
1.4xl0 -2 
0.13 

6.1xl0 -2 

3.3xl0~ 2 

0.09 

26.6xl0 -2 

1.7xl0 -2 

0.04 


The methods applied for the removal of the different 
perturbation sources are briefly discussed in the following 
section. 

8.1.4 Estimate of Measurement Perturbations 

The tower-shading correction factors ^^(A), (A), 

and rf^ s (X) were computed using a Monte Carlo modeling 
of the radiance and irradiance fields at the AAOT site 
for the specific point and geometry where the WiSPER 
measurements were taken. Modeling results were validated 
through field data obtained from an experiment specifically 
designed to quantify L u (z, A) and P<f(z,A) tower-shading 
effects as a function of the distance from the tower (Zibordi 
et ah 1999). 

The operational tower-shading corrections (Doyle and 
Zibordi 2002) were performed making use of look-up ta- 
bles (LUTs) in which the correction factors for each center 
wavelength are indexed by discrete values of: 

■ The solar zenith and azimuth angles, 9q and 0o, 
respectively; 

■ The total seawater absorption, a( A); 

■ The seawater single scattering albedo, defined as 
lj 0 (A) = b(A)/[a(A) + b( A)], where b( A) is the total 
seawater scattering coeffficient; and 

■ The above- water diffuse over direct irradiance ratio, 
r z ( A) . 

In building the LUTs, the variability in r*( A) was produced 
by simulating the tower-shading uncertainties during clear- 
sky conditions with different values of the aerosol optical 
thickness, r a ( A). The overcast sky condition was simulated 
as a special case by assuming a totally isotropic sky illu- 
mination, i.e., equivalent to the light field resulting from 

To( A) -> 00. 

Table 3 presents the discrete reference values 9 0 , 0o, 
a( A), cjq(A), and r a (A) used for computing the LUTs for 


the tower-shading correction scheme. The OCI-200 and 
OCR-200 wavelength-dependent values for the bottom re- 
flectance pb( A) 7 Rayleigh (molecular scattering) optical 
thickness tr( A), and ozone optical thickness t 0 { A), are 
given in Table 4. The discrete reference values shown in 
Table 3 were chosen based on the representative variability 
of the related quantities measured at the AAOT site. 

The a (A) values used for retrieving correction factors 
from the LUTs are computed from AC-9 measurements 
or, as an alternative when AC-9 data are not available, ob- 
tained from discrete surface seawater samples. Similarly, 
the b( A) values are computed from AC~9 measurements or, 
as an alternative, estimated from the subsurface irradiance 
reflectance, R((T, A) = £ u (0~, A)/&(<r, A), using the re- 
lationship suggested by Kirk (1994): 


b( A) = o( A) 


103P(Q~,A) 

1 - R{ 0",A)’ 


(7) 


where the tower-shading, self-shading, and bottom effects 
on R( 0“,A) are neglected. 

The possible sources of uncertainties in the correction 
scheme for removing tower-shading effects, some of which 
might be significant under certain circumstances, are as 
follows: 

1) Assuming the tower structure absorbs all incident 
photons (i.e., the tower is perfectly black and does 
not reflect any light); 

2) Neglecting roughness in modeling the sea surface 
reflectance; and 

3) Using correction factors computed for discrete val- 
ues of the input variables 0o, (f> o, a(A), u> 0 (A), and 
r a (A) without any interpolation. 

The instrument self-shading correction factors *7^ (A) 
and tT? 5 ( A) are computed using the scheme proposed by 
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Gordon and Ding (1992) plus the parameterizations sug- 
gested by Zibordi and Ferrari (1995) and Mueller and Aus- 
tin (1995). Setting the OCI-200 and OCR~200 instrument 
radius Rj = 4.5 cm and the radius of the sensor entrance 
optics R s = 0.5cm, the computation of the self-shading 
correction factors is made using the following data: 

a) The average i9 0 of each cast; 

b) a( A) computed from AC-9 measurements or as an 
alternative (when AC-9 data are not available) ob- 
tained from surface seawater samples; and 

c) r t ( A) from the MFR-6 data taken at the same time 
as the in-water optical profiles. 

Possible sources of uncertainties in the computation of 
7 lfS(X) and 7}£ S (A) include the following: 

■ Assuming the instrument as an ideal disc instead of 
a cylinder; 

■ Assuming each channel sensor is located at the cen- 
ter of the instrument (they are actually arranged in 
a circle with one channel in the center); and 

■ Neglecting the roughness of the sea surface. 

The correction factors r?|^(A) and vf E W for bottom 
(sand and mud at approximately 17m depth), are com- 
puted making use of a simple analytical model derived from 
Maritorena et al. (1994) for bottom correction of irradiance 
reflectance. The model requires the so-called operational 
diffuse attenuation coefficient, KE d ( A), which is assumed 
here to be equal to the average diffuse attenuation coef- 
ficient obtained from the linear regression of In \Ed{z, A)] 
versus depth between the surface and the bottom depth, 
zq. Another needed parameter is the bottom reflectance, 
Pb{ A), which is given by ( 


Similar to E EE ( 0 ,A), the bottom-corrected upward 
radiance L^((T,A) is estimated by: 


L* E ( O’, A) = 4(0“, A) 


Pb(A)4(0 ,A)k(z B ,A) 

Qn( 0-,A) 


where, as given in (1), Q n (0“,A) = £ U (<T, A)/4(0", A) 
in units of steradians. The former relationship assumes 
Qn(0',A) is not affected by bottom effects, i.e., E u ( O', A) 
and L u ( 0", A) are identically affected by the bottom. Even 
though this assumption is only valid in isotropic conditions 
[i.e., when Q n (0 , A) = 7rsr and the bottom reflectance 
Pb( A) is lambertian] it can be shown that in the worst 
conditions observed in the entire CoASTS data set (i.e., 
with vf E W « 0-9), an extreme uncertainty of ±10% in 

Qn(0“, A) induces a variation of approximately ±0 5% in 

£?*(0-,A). 

Uncertainties in the computation of rj£ E ( A) and rj? E ( A) 
are induced by: 

a) Applying a simple (approximate) analytical model; 

b) Using the operational diffuse attenuation coefficient- 
and 

c) Assuming the bottom reflectance is lambertian. 

8.1.5 Environmental Effects 


Pb{ A) 


&u(zB) A) 
4(^b, A) 


where E u (z B , A) and Ed(z B , A) are the extrapolated values 
E u (z, A) and Ed(z i A) at the bottom depth z B (the mean 
Pb{ A) values at the AAOT site are listed in Table 4). 

The bottom-corrected upward irradiance just below the 
surface is given by 


? (0~, A) = 


-k{z B , a) r u(0 > A) 

Pb(A) E d (0-, A) k(z B , A)j , 


where k(z Bl A) is the transmittance for the downward plus 
upward normal optical paths beween the surface and the 
bottom at depth z B and is equivalent to 

k ( z B, A) = exp [—2 K Bd (A) z B ] . (10) 


To complete the analysis on radiometric quantities, the 
uncertainties due to surface roughness, changes in sea- 
water characteristics, and variability in the illumination 
conditions during the in-water optical profiling are inves- 
tigated. These uncertainties are estimated for L u (0',A) 
and E d (0 , A) using data taken during clear-sun conditions, 
with wind speeds less than 5 ms -1 (the typical circum- 
stance for CoASTS measurements), and at a time close to 
local noon (to minimize variations in illumination condi- 
tions caused by changes in the solar zenith angle). 

Using 19 stations that fulfill the aforementioned condi- 
tions, the uncertainty attributed to environmental effects 
was estimated from the average absolute percent difference) 
(APD) between L u { O’, A) values, and between E d ( O', A) 
values, computed from successive profiles (i.e., the first 
and the third profile of each measurement station) about 
20 min apart. Results for L„(0',A) show average APDs 
ranging from 2.7 ±3.1% to 3.2 ±2.9% at 412 and 665 nm, 
respectively. Similarly, results for £^(0",A) show average 
APDs ranging from 1.9 ± 1.7% to 3.0 ±2.3% at the same 
wavelengths. The former relatively low values are probably 
the consequence of two factors: 

f The APD between two quantities, X and Y, is defined as 
200 1 A - Y\/(X + Y). 
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Table 5. Average values of r t (A) measured at the AAOT site during clear-sun conditions. The standard 
deviations are indicated by the “±” sign. — . 1 


aeviauunb aie 

A [nm] 

412 

V, _1_ 

443 

490 

510 

555 

665 

Average 

0.89 ±0.41 

0.73 ±0.35 

0.53 ±0.27 

0.46 ±0.24 

0.34 ±0.19 

0.23 ± 0.14 


a) The high stability of the deployment system, which 
prevents any sensor tilts; and 

b) The low winch deployment speed (0.1 m s 1 ), which 
ensures at least 60 samples per meter for the aver- 
aging of surface perturbations in the profile data. 

8.2 Ed( 0 + ,A) and £;( 0%A) Measurements 

Total and diffuse above-water solar irradiances are mea- 
sured using the MFRr6 automatic radiometer (Harrison et 
al. 1994). Measurements are performed in six bands 10 nm 
wide (Table 2) and in one broad band (ranging from ap- 
proximately 400— 1 ,100 nm) . The total and diffuse compo- 
nents are measured by alternately exposing and shading 
the cosine collector of the instrument. The collector shad- 
ing from the direct sun component is obtained through an 
automated arc-shaped (shadow) band moving above the 
center of the entrance aperture, which blocks a portion of 
the sky with a 3.3° angle. 

Each measurement sequence includes four independent 
measurements: 

■ Ed{ 0 + ,A) taken with the band at its home (i.e., 
nadir) position; 

■ Ei( 0 + , A) with the sun blocked by the band; 

■ Two more measurements taken when the band is 
positioned 9° to either side of the sun. 

The latter two measurements are automatically used by 
the system firmware to compute Ei( 0 + ,A) by correcting 
Ei{ 0 + ,A) for the sky radiance blocked by the band when 
the sun and a portion of the sky are occulted. The system 
automatically performs a maximum of three sequences of 
measurements per minute. 

In the framework of the CoASTS campaigns, MFR-6 
data were taken at 10 min intervals during each measure- 
ment station. The data, used to estimate 7q(A) for self- 
shading and tower-shading corrections, are computed ac- 
cording to 


r,W E d { 0 + ,A) - Ei( 0 + ,A)' 


(12) 


8.3 E s ( A) and X) Measurements 

The direct normal sun irradiance, E S (X) and the diffuse 
sky radiance, in a wide range of angles in the almucantar 
and sun planes, are both measured with a CE-318 auto- 
matic sun photometer. The device is composed of. 

1) A sensor installed in an alto-azimuthal platform, 

2) A programmable unit controlling the measurement 
sequences, sun— sky pointing, and data logging, and 

3) A data transmission unit based on the Meteorolog- 
ical Satellite (METEOSAT) Data Collection Plat- 
form (DCP) system. 

The optical part of the CE-318 is composed of two col- 
limators with a 1.2° full-angle field of view [one is used 
for Li(9,cj>,X) and the other is used for both E( A) and 
Li{0, 4>, A) in the sun aureole], and a filter wheel with eight 
filters. The direct normal sun irradiance, used for the re- 
trieval of the aerosol optical thickness, is taken at regu- 
lar air mass intervals in all spectral channels. The diffuse 
sky radiance, mainly used for the retrieval of the aerosol 
particle size distribution and phase function (Holben et 
al. 1998), is only taken for channels in the 440-1,020 nm 
range (Table 2). The CE-318 installed on the AAOT, is 
part of the Aerosol Robotic Network (AERONET) at the 
NASA Goddard Space Flight Center (GSFC) where the 
regular instrument calibrations occur , as well as the almost 
near-real time data processing and data quality assurance 
(Holben et al. 1998). 

Instrument calibrations for Li(9,4>, A) data were per- 
formed by the AERONET project using a 2 m integrating 
sphere. The uncertainty in the absolute radiometric cali- 
bration of L t {6 , 4>, A) is assumed to be within ±4%. Instru- 
ment calibrations for E( A) data were carried out through 
intercalibration with a reference sun photometer regularly 
calibrated at the high altitude Mauna Loa Observatory 
(Hilo, Hawaii) with the Langley technique (Holben et al. 
1998). 

The calibrated sun photometer data E( A), assuming 
no water vapor absorption, are related to the atmospheric 
optical thickness, r, through 

E( A) = Eo{X)D y exp[— t(A) m], (13) 


Because E d { 0 + ,A) and £,(0 + , A) are measured with the 
same instrument, an absolute radiometric calibration of 
the instrument is not required. Table 5 gives the values 
for ri( A) measured at the same time as the in-water op- 
tical data during clear-sun conditions for the 1995-1998 
period. Under fully cloudy conditions, r t (A) -> oo. 


where E 0 { A) is the extra atmospheric irradiance, D y is the 
sequential day of the year (which accounts for the sun- 
Earth distance), r = T/j+T 0 +r a (i.e., the sum of the molec- 
ular, ozone, and aerosol optical thicknesses, respectively), 
and m is the relative air mass. Assuming the instrument 
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is properly calibrated, the uncertainty associated with the 
retrieved r a is estimated to be less than 0.02 (Smirnov et 
al. 2000). 

Results from Langley calibrations of E( A) data taken at 
the AAOT on clear and atmospherically stable days, have 
shown differences generally lower than ±2% with respect to 
the calibrations performed by GSFC in the 440-1,020 nm 
spectral range. This result supports the concept of using 
clear and stable days for making on-site field checks of the 
stability of the absolute radiometric calibration. 

8.4 c{z, A) and a(z,\) Measurements 

The beam attenuation c(z, A) and absorption a(z, A) 
coefficients of seawater (both in units of per meter) are 
computed from measurements taken in nine bands (Ta- 
ble 2), 10 nm wide, with a 25 cm pathlength AC-9 equipped 
with a Sea-Bird Electronics (Bellevue, Washington) 5T 
submersible pump for constant water flow. The calibrated 
beam attenuation and absorption coefficients, c t „ w (z,\) 
and a t - w (z, A), respectively, are obtained from the AC-9 
measurements for suspended and dissolved optical compo- 
nents (not including the contribution of pure seawater). 
Both of these parameters are then corrected for tempera- 
ture, T w , and salinity, S w , differences between the in situ 
seawater and the pure water used for laboratory calibra- 
tion (following the methodology recommended in an inter- 
nal document by WETLabs, the instrument manufacturer, 
with CTD profile data taken at the same time as the AC-9 
measurements). 

While c?T w {z,\) — the measured beam attenuation co- 
efficient corrected for salinity and temperature effects — 
does not require any further processing [i.e., c(z, A) = 
cfj w (z 1 A) -bc^A), where ^(A) is the pure seawater beam 
attenuation coefficient] , the resulting salinity and tempera- 
ture corrected afl w {z } A) data need to be further corrected 
for scattering effects to retrieve a(z, A). The latter scatter- 
ing effects, due to the finite acceptance angle of the optics 
and to the incompletely reflective surface of the absorption 
chamber (which prevents the detector from collecting all 
the scattered light), induces overestimates of the retrieved 
absorption coefficient. 

The WETLabs suggestions for correcting the absorp- 
tion data for scattering effects, recounted here for com- 
pleteness, are based on three methods: 

1- The removal at each wavelength of the measured 
absorption coefficient at a reference wavelength, A 0 , 
in the near infrared part of the spectrum (usually 
Ao = 715 nm) where the absorption by particulate 
matter and colored dissolved organic matter is as- 
sumed to be zero: 

a(*,A) = afI w (z,X) + a w (\) - afl w (z, A 0 ), (14) 

where a^(A) (in units of per meter) is the pure sea- 
water absorption (Pope and Pry 1997). 


2. The removal of a fixed percentage, 5 (which is inde- 
pendent of A and typically varying from 0.14-0.18), 
of the scattering coefficient 6 t _ w (z, A) = c t _ w (z, A)- 
a t-w(z, A), by assuming that the shape of the vol- 
ume scattering function is independent of A and of 
the type of particulate material: 

o(z,A) = afZ w (z, A) + a w (X) 

~S[cfT w (zA)-&fI w (zA)l (15) 

3. The removal of a variable percentage of the scatter- 
ing coefficient b t - w (z, A) estimated from cfJ^(z,A) 
and afjf^^A) at A and A 0 , by assuming that the 
absorption coefficient of particulate and dissolved 
material is zero at the reference wavelength Aq — 
715 nm and that the shape of the volume scattering 
function is independent of wavelength (Zaneveld et 
al. 1994): 

a(z,A)= af^ w (z, \) + a w (\) 

(l 6 ) 

_ Aq) [cfjjz, A) - afjjz, A)] 

cfl w (z, A 0 ) - afj w (z, A 0 ) 

In Fig. 4, an intercomparison is given between the sum 
of the absorption coefficients obtained from discrete sea- 
water samples (Sects. 8.6 and 8.7), a ph {z , A) + a dp (z , A) + 
a ys( z . A), and the absorption coefficients computed from 
the AC-9 data corrected for scattering effects according to 
method 3 above, a t - w (z, A) = a(z, X)-a w {z, A). The inter- 
comparisons are given for the 412, 488, 555, and 676 nm 
wavelengths, at depths z= 1, 8, and 14 m for which ab- 
sorption coefficients from the AC-9 and seawater discrete 
samples are both available. The AC-9 values are for data 
binned with aim depth resolution. The agreement be- 
tween a ph (z, A) + a dp (z, A) + a y9 (z, A) and a t - w (z. A), is 
very good at 412 nm, but a general underestimate of the 
absorption is noticed at 488, 555, and 676 nm with the 
worst results at 510 nm (not shown here), which exhibit 
average percent differences of 25%. 

The absolute calibration of the AC-9 was performed 
by WETLabs. Additional calibrations were periodically 
performed in the laboratory using clean (Milli-Qf) wa- 
ter. Offsets between the original WETLabs and labora- 
tory calibrations were used for correcting c t _ w (z,\) and 
at-w(z> A) measurements by drifts due to changes in in- 
strument sensitivity. The former clean-water calibrations 
showed offsets slightly drifting over the measurement pe- 
riod and exhibiting values generally from 0.08-0.14 m" 1 for 
c t ~ w (z, A), and from 0.005-0.07 m _1 for a t ~ w (z, A) in the 
spectral range 412-715 nm. Recent work on AC-9 calibra- 
tion (Twardowski et al. 1999), has shown the need for daily 
clean-water calibrations. Because of this, it may result that 

t Milli-Q is a trademark of Millipore Corporation (Bedford, 
Massachusetts) . 
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o t _ w [rrT 1 ] °t-w [m ] 


Fig. 4. Scatterplots of a ph + a dp + a ys from discrete seawater samples and a t - w from AOS I data computed 
by applying the scattering correction method 3 (using 177 samples at 1, 8, and 1 m). e as e lne 
(whose intercept A, slope B, and the determination coefficient r 2 are given in the panels show the fitting 
line produced with the “robust” least absolute deviation method (Press et al. 1987). Plots are given for 
wavelengths of a) 412 nm, b) 488 nm, c) 555 nm, and d) 676 nm. 


the clean-water calibrations performed in CoASTS during 
the 1995-1998 period, do not fully resolve short-term drifts 
in the instrument response. 

8.5 Ancillary Field Data Measurements 

Salinity S„,(z) (in practical salinity units) and temper- 
ature T w (z) (in units of degrees Celsius) profiles are ob- 
tained at the same time as the WiSPER and AC-9 profiles 
(and are needed for AC-9 data processing). The tide level 
Ti and the Secchi depth SD were recorded for each sta- 
tion, together with major meteorological quantities taken 
at about 15 m above the sea surface: atmospheric pressure, 
P a (in units of hectopascals); relative humidity, RH (in 
percent); air temperature, T a (in degrees Celsius), wind 


speed, W a (in meters per second); and wind direction, 
Wd (in true degrees). In addition, sea state, M [using 
the World Meteorological Organization (1983) code], and 
cloud cover, CC } observations in quarters are recorded 
to provide qualitative information for interpreting (appar- 
ently) spurious cases. 

8.6 a p (A) Measurements 

The in vivo absorption coefficient, a p (A), of aquatic 
particles retained on filters (in units per meter), is deter- 
mined with a Perkin Elmer (Fremont, California) Lambda- 
19 dual-beam spectrometer equipped with a 60 mm diam- 
eter integrating sphere. The deposit of the particles on fil- 
ters is obtained by filtration of seawater samples on glass 
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fiber (GF) filters with a nominal pore size of 0.7 ^m under 
low vacuum pressure (less than 120mmHg) to prevent par- 
ticle breakage and pigment degradation. The filtered vol- 
ume varies from 1.0-2. 5 L, and is a function of the quantity 
of material suspended in the seawater. After filtration, the 
filters with the deposits are immediately placed on a Petri 
slide and stored in liquid nitrogen. The total absorption 
coefficient of the equivalent particle suspension in the 400- 
750 nm spectral range (with 1 nm resolution) is computed 
according to: 

o p (A) = 2.3A SUS (X)^, (17) 

where V w is the volume of filtered water in liters, F a is the 
filter clearance area in square centimeters, and A JUS (A) 
is the equivalent particle suspension absorbance obtained 
from the Transmission and Reflection (T-R) method pro- 
posed by Tassan and Ferrari (1995), which has been shown 
to be appropriate for the analysis of water samples char- 
acterized by highly backscattering mineral particles or by 
highly absorbing sediments. 

Alternative to the T-R method, is the transmission 
method (T), which simply assumes that particle absorption 
is negligible in the near infrared and the effects of scatter- 
ing are wavelength independent (Mitchell and Kiefer 1988). 
These simpler assumptions may be valid for oceanic wa- 
ter dominated by phytoplankton exhibiting no absorption 
in the near infrared, but are inadequate for coastal water 
containing large fractions of nonpigmented material (which 
may show an appreciable absorption in the near infrared, 
as well as wavelength-dependent scattering). 

The two components of the particulate absorption co- 
efficient for the pigmented and nonpigmented fractions, 
a p h(A) and dd p ( A), respectively, are obtained through 
bleaching the sample on the filter using a solution of sodium 
hypochlorite (NaCIO) as an oxidizing agent. This oxida- 
tion acts rapidly on pigment molecules and slowly on detri- 
tus making possible a selective analysis of the absorption 
components of pigmented and nonpigmented particles re- 
tained on the filter. A description of the NaCIO bleaching 
technique, is presented in Tassan and Ferrari (1995) and 
in Ferrari and Tassan (1999). 

Uncertainty in the estimate of the particle absorption 
coefficients is induced by the use of GF filters with a nom- 
inal 0.7/im pore size. In fact, these filters do not allow 
bacteria and the fraction of mineral particles with diam- 
eters lower than 0.7^m to be accounted for. In general, 
however, the absorption of these small mineral particles 
is negligible compared to the total absorption. The ab- 
sorption of bacteria is almost 10 times lower than that of 
algal cells, and 5-10 times lower than that of cyanobac- 
teria (Morel and Ahn 1990) characterized by a mean cell 
diameter of approximately 1 /x m. 

Figure 5 shows the scatterplots of a p (A) obtained with 
the T method versus the T-R method at the representa- 
tive wavelengths 412, 490, 555, and 665 nm. The a p (A) 


data in Fig. 5 shows a good agreement at 555 nm with 
average percent differences of 5% displayed by the least- 
squares fit. The a p (A) values obtained with the T method, 
however, show some overestimate with respect to the T-R 
method at wavelengths shorter than 510 nm (exhibiting av- 
erage percent differences up to 11%) and some underes- 
timate at longer wavelengths (exhibiting average percent 
differences up to —15% at 555 nm). The latter differences 
are due to the assumptions of a p (750) = 0 and wavelength- 
independent of scattering in the T method). Pronounced 
differences can be observed (not shown here) between the 
a p h( A) data obtained with the two methods applied to 
the bleached samples, while the a p / l (A) data from the T 
method only exhibit a slight general overestimate with re- 
spect to the T-R method. The dependence on wavelength 
of the regression fits given in the different scatterplots in 
Fig. 5 is a clear indication of the relevance of the spectral 
dependence of the scattering correction for the considered 
samples. 

The repeatability and reproducibility of the CoASTS 
in vivo particulate absorption measurements were investi- 
gated using two types of analyses: 

1) Duplicate analysis of the same samples (each filter 
sample was analyzed twice) for repeatability; and 

2) Analysis of replicate samples (two separate samples 
obtained from the same water collection event were 
analyzed) for reproducibility. 

Results from the two types of analyses are summarized in 
Table 6 at different wavelengths. 

The duplicate analysis of the same samples, show aver- 
age APDs varying from 3.0 ±2.6% at 412 nm for a p (A) = 
0.084 ± 0.044 m” 1 , up to 7.4 ± 6.0% at 555 nm for a p (A) = 
0.023 ±0.011 m -1 . The former differences are attributed to 
instrument sensitivity, which can induce large uncertain- 
ties in the presence of samples characterized by weakly 
absorbing deposits on filters (e.g., because of the filtration 
of a small seawater volume). A second reason for the for- 
mer differences could be a slight variation in the mechan- 
ical repositioning of the sample in front of the aperture 
of the integrating sphere, added to the spatial nonhomo- 
geneity of the deposit on the filter. These combined effects 
may significantly change the sample area viewed by the 
spectrometer through the integrating sphere aperture of 
approximately 150 mm 2 . 

The analysis of replicate samples shows average APDs 
of 9.5 ± 5.5% at 412 nm for a p (A) — 0.093 ± 0.051m -1 , and 
9.8 ±7.0% at 555 nm for a p (A) = 0.024 ± O.O^m^ 1 . The 
former differences, increased by a few percent with respect 
to those given for the duplicate analysis of samples, are 
justified by 

■ The spatial nonhomogeneity of the deposits on the 
filters, 

■ Actual differences between samples, and 

■ Instrument sensitivity. 
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Fig. 5. Scatterplots of a p obtained with the T method versus the T-R method (using 332 surface samples). 
The dashed line indicates the least-squares fitting line (whose intercept A, slope B, and the determination 
coefficient r 2 , are given in the panels). Plots are given for wavelengths of a) 412 nm, b) 490 nm, c) 555 nm, 
and d) 665 nm. 


Table 6. Summary results on the analysis of a p (A) measurements at different wavelengths (in nanometers). e 
average APDs resulting from duplicate analysis of the same sample is d\ . The average APDs resulting from anajsis 
of replicate samples is d 2 - The parentheses indicate the average a p (A) values (in units of per meter) for the analyzed 
samples, and the standard deviations are indicated by the “±” sign. 


3.0 ± 2.6 2.9 ±2.3 3.7±2.7 

(0.084 ±0.044) (0.082 ±0.042) (0.055 ±0.026) 


8.9±5.9 8.3 ± 5.4 

(0.090 ±0.049) (0.059 ± 0.030) 


4.6 ±3.3 7.4 ±6.0 

(0.041 ±0.019) (0.023 ±0.011) 


8.3 ±5.2 
(0.044 ±0.021) 


9.8 ±7.0 


7.2 ±5.8 
(0.027 ±0.015) 

10.7 ±7.6 


(0.024 ±0.012) I (0.028 ±0.017) 


L5 
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Table 7. A summary of the results on the analysis of a ys ( A) measurements at different wavelengths (in 
nanometers). The average APDs resulting from duplicate analysis of the same sample is d\, and the average 
APDs resulting from the analysis of replicate samples (in this case, two replicates) is d 2 . The parentheses 
indicate the average a yfl (A) values (in units of per meter) for the analyzed samples, and the standard deviations 
are indicated by the “±” sign. 


A [nm] 

412 

443 

490 

510 

555 

di [%} 
n = 21 

10.1 ±7.3 
(0.168 ±0.037) 

12.4 ±8.4 
(0.097 ±0.022) 

20.7 ±15.0 
(0.044 ±0.011) 

21. 5± 15.7 
(0.031 ± 0.007) 

24.2 ±19.8 
(0.015 ±0.005) 

d 2 [%] 

n = 21 

12.1 ±6.3 
(0.175 ±0.038) 

15.8 ±8.8 
(0.103 ±0.026) 

26.9 ±17.5 
(0.047 ±0.015) 

25.8 ±18.0 
(0.035 ±0.011) 

30.3 ±23.8 
(0.018 ±0.005) 


The study on reproducibility of a p (A) also gives an esti- 
mate of the precision of the whole measurement procedure 
applied in CoASTS (including seawater sampling, filtration 
of samples, and spectrometric analysis). Focused studies 
on the accuracy of the T-R method are given in Tassan 
and Ferrari (1995 and 2002). 

8*7 a ys ( A) Measurements 

The absorption coefficient of CDOM or so-called yel- 
low substance , a y5 (A) in units of per meter, is determined 
using a dual-beam spectrometer Perkin Elmer Lambda- 
12. Water samples (generally 0.7-1 L) are filtered through 
0.22 /xm cellulose filters and then refrigerated at 4°C in an 
amber glass bottle with the addition of a solution (0.5 mL 
for 100 mL of sample) of 10 gL” 1 of sodium azide (NaNa) 
to prevent the growth of any bacteria. The spectrometric 
measurements, which are generally carried out within a few 
days, are performed in the spectral range 350-750 nm with 
1 nm resolution. They are performed by placing a 10 cm 
quartz cuvette containing Milli-Q water in the optical path 
of the reference beam, and a 10 cm quartz cuvette contain- 
ing the filtered seawater sample in the optical path of the 
sample beam. The spectral absorption coefficient a ys ( A) is 
computed from the measured absorbance A y3 (A) resulting 
from the difference between the sample absorbance and the 
reference absorbance (Ferrari et al. 1996): 

dy 3 (X) = 2.3 ( 18 ) 

where A y3 (A) is the measured absorbance, and L c is the 
pathlength (in meters) of the cuvette. The instrument 
background is removed using measurements obtained with 
Milli-Q water in both the sample and the reference cu- 
vettes. 

The use of 0.22 /xm pore size filters to define CDOM, 
when 0.7 /xm pore size filters are used for particle absorp- 
tion, indicates the overall absorption budget cannot be 
fully resolved. In fact, bacteria, colloids, and very small 
mineral particles having sizes between 0. 2-0.7 /xm are not 
included in the absorption analysis. Their inclusion in 
the filtered seawater sample for a ys (A) analysis, however, 


would likely lead to an overestimate of the absorption co- 
efficient because of their very high scattering and low ab- 
sorbing properties (Stramsky and Mobley 1997, and Pak 
et al. 1971). 

The repeatability and reproducibility of a ys (A) mea- 
surements were investigated through duplicate analysis of 
the same samples and analysis of replicate samples, re- 
spectively. The analysis was carried out using a ys ( A) spec- 
tra corrected by the background signal at 600 nm [i.e., for 
each spectrum, a ys (600) is subtracted from a ys { A) at all 
wavelengths assuming CDOM is not absorbing in the red 
domain]. The results of this analysis are summarized in 
Table 7 at different wavelengths. 

The duplicate analysis of the same samples shows aver- 
age APDs vary as a function of the absorption value from 
10.1 ± 7.3% at 412nm for a ya (A) = 0.168 ± 0.037m -1 , up 
to 24.2 ± 19.8% at 555 nm for a ya (A) = 0.015 ± 0.005 m” 1 . 
These differences are mostly attributed to temperature dif- 
ferences between the seawater sample and the reference 
Milli-Q water, the formation of microbubbles on the wall 
of the cuvettes, and mechanical alignment of the cuvettes. 

The analysis of the replicate samples shows increased 
average APDs, when compared to the duplicate analysis, 
varying from 12.1 ±6.3% at 412nm for a y5 (A) =0.175 ± 
0.038 m“ 1 up to 30.3 ± 23.8% at 555 nm for a y5 (A) = 0.018± 
0.005 m” 1 . The latter increase in differences, with respect 
to differences computed from duplicate analysis, are justi- 
fied by differences between samples. 

8.8 Pigment Measurements 

Phytoplankton pigment concentrations are determined 
using HPLC analysis with a slightly modified Joint Global 
Ocean Flux Study protocol ( JGOFS 1994 and Jeffrey et al. 
1997). Seawater volumes, ranging from 1.0-2. 5 L, set as a 
function of the quantity of material suspended in the sea- 
water, are filtered immediately after collection using GF 
filters with a nominal pore size of 0.7 /xm. The filters are 
then immediately stored in liquid nitrogen for subsequent 
laboratory analysis. Using the 0.7 /xm pore size is justi- 
fied by the diameter of living phytoplankton cells, which 
is generally higher than 1 /xm (Stramsky and Kiefer 1991). 

Pigment extraction is carried out in dim light by soak- 
ing each filter in 5mL of a solution composed of 100% 
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acetone and 0.2 /xgL -1 canthaxanthin [used as an internal 
standard (IS)], and then grinding it with a Teflon® f pestle 
for 30 s at 0°C in a thermostatic tube. After storing the 
reworked sample for 24hina— 18°C freezer, it is filtered in 
dim light using a syringe with a 0.45 /xm GF filter and then 
taken to ambient temperature. Finally, 150 /xL of a buffer 
solution composed of 0.5 M ammonium acetate (which en- 
sures a better peak separation) is added to 500 /xL of the 
filtrated extract and injected into the 200 /xL HPLC loop. 
From the HPLC chromatogram, the concentration of each 
pigment, C p (/xgL -1 ), is computed using 


The scatterplot of total pigment concentration versus 
total chlorophyll a (chlorophyll a plus chlorophyllide a), is 
shown in Fig. 6 as an empirical quality assurance of the 
data. The results exhibit a high correlation (98% of the 
variance explained by the linear fit), and the slope of the 
linear fit (i.e., 2.05) is in close agreement with the value 
of 1.96 published by Aiken et al. (1995) for a data set 
including measurements from different oceanographic bio- 
optical provinces. 


C p — 


A p f p Ve 

V t V w Be 


Ar 

A. 


(19) 


where A v [in milli- Absorbance Units (AU) per second] is 
the pigment peak area; f p is the relative response factor 
for each pigment; V e (in milliliters) is the volume of the 
pigment extract; Vi (in microliters) is the volume injected 
into the HPLC; V w (in liters) is the filtered volume; B e is 
the buffer extract dilution factor (i.e., the ratio between 
the 500 fjJu of extract volume and the sum of the 500 /xL 
of extract volume to the 150 /xL of buffer solution); A r (in 
milli- AU per second) is the reference IS peak area; and A s 
(in milli-AU per second) is the sample IS peak area. The 
factor f p (in units of ngmAU^s -1 ) is specific for each 
pigment and is obtained through calibration of the instru- 
ment response with pigment standards. The reference IS 
peak area, A r , is computed from the HPLC chromatogram 
of a reference solution of 0.2 ^xgL” 1 IS ( A r is recomputed 
at least once per day). 

The use of an IS allows correction for uncertainties 
caused by evaporation or experimental losses during pig- 
ment extraction. System calibration is performed on an 
annual basis using pigment standards provided by the DHI 
Water and Environment Institute! (H0sholm, Denmark). 
The list of analyzed pigments are as follows: 

■ Chlorophyll a, 

■ Chlorophyll 6, 

■ Chlorophylls ci and C2, 

■ Chlorophyllide a, 

■ Fucoxanthin, 

■ Diadinoxanthin, 

■ /3-carotene, 

■ Zeaxanthin, 

■ Alloxanthin, 

■ 19'-butanoyloxyfucoxanthin, 

■ lO'-hexanoyloxyfucoxanthin, and 

■ Diatoxanthin. 



f Teflon is a registered trademark of E.I. du Pont de Nemours, 
Wilmington, Delaware. 

\ Formerly the VKI Water Quality Institute. 


Fig. 6. A scatterplot of total pigments versus to- 
tal chlorophyll a ( C a )> defined as chlorophyll a plus 
chlorophyllide a (including intercept A, slope B, 
and the determination coefficient r 2 ). 

In order to evaluate the precision of measurements per- 
formed with the described method, different experiments 
were carried out to estimate repeatability through dupli- 
cate analysis of the same samples and reproducibility using 
two methods: 

a) Analysis of replicate samples, and 

b) An intercomparison analysis of replicate samples by 
a different laboratory using a different measurement 
system and a different method. 

Results from these experiments are summarized in Table 8 
for some of the most represented pigments in the northern 
Adriatic Sea. 

The duplicate HPLC analysis of samples shows aver- 
age APDs in C a of 1.4 ± 1.2% for 24 samples with an av- 
erage concentration of 1.18 ± 0.65 /xgL" 1 (this variability 
can only be attributed to the sensitivity of the measure- 
ment system). The analysis of replicate samples, shows 
average absolute percentage differences in C a of 7.1 ± 6.3% 
estimated using 36 pairs of samples with a mean concentra- 
tion of 0.86 ±0.19 /xgL -1 (these differences are attributed 
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Table 8. A summary of the results on the analysis of HPLC data for the most represented pigments in the 
northern Adriatic Sea. The average APDs resulting from the duplicate analysis of the same sample is d\. The 
average APDs resulting from the analysis of replicate samples is d 2 . The average APDs from replicate sample 
analyses performed by JRC and by LOV is The values in parentheses indicate the average concentration of 
pigments (in units of microgram per liter) for the analyzed samples, and the standard deviations are indicated 
by the “±” sign. 


Pigment 

Chlorophyll a 

Chlorophyll b 

Chlorophyll c\ + c 2 

Fucoxanthin 

1 9 '-Hexanoyloxy- 
fucoxanthin 

di [%] 
n = 24 

1.4 ± 1.2 
(1.18 ±0.65) 

2.8 ±3.0 
(0.07 ±0.05) 

4.7 ±5.7 
(0.12±0.10) 

1.1 ±1.5 
(0.42 ±0.30) 

3.2±3.1 
(0.14 ±0.11) 

d 2 [%] 
n = 336 

7.1 ±6.3 
(0.86 ±0.19) 

7.8 ±6.9 
(0.13±0.10) 

14.4 ±10.8 
(0.08 ±0.03) 

5.2 ±4.1 
(0.14±0.12) 

8.1 ± 7.1 
(0.08 ±0.04) 

d 3 [%] 

n = 15 

13.1 ±7.0 
(2.05 ±0.97) 

14.9 ±16.7 
(0.12 ±0.07) 

30.2 ±8.8 
(0.49 ±0.33) 

21.0 ±8.8 
(0.98 ±0.63) 

5.4 ±13.1 
(0.17 ± 0.13) 


to the sensitivity of the measurement system, to the ef- 
ficiency of the pigment extraction, and to differences be- 
tween the duplicate samples). 

The intercomparison experiment, which used sample 
analyses from a different laboratory, involved the LOV, 
and was carried out in September 1998. Pigments were 
extracted by applying the methods previously described. 
The analysis was performed with different systems and 
by applying different methods to replicate samples. The 
CoASTS (JRC) method was the same one as described 
above, and LOV used the method described in Vidussi et 
al. (1996). The results showed average APDs in C a of 
13.1 ±7.0% estimated over 15 samples with a mean con- 
centration of 2.05 ± O^pgL -1 . These differences are at- 
tributed to the use of different HPLC methods, different 
sensitivities of the measurement systems, different calibra- 
tion accuracies, and differences between the duplicate sam- 
ples. 

An estimate of reproducibility between different labo- 
ratories is given by the difference between the uncertainty 
associated with the analysis of the replicate samples (d 2 ) 
and the total uncertainty estimated from the laboratory in- 
tercomparison (d 3 ). This reproducibility estimate (about 
6% for chlorophyll a) is very similar to the reproducibility 
achieved during a recent HPLC intercomparison (Hooker 
et al. 2000b) based on open ocean samples (on average 
about 5.5% for chlorophyll a). 

8.9 CpsM Measurements 

TSM concentration is obtained from the net weight of 
the material collected on GF filters with a 0.7 pm pore 
size following a slightly modified version of the method 
proposed by Strickland and Parsons (1972). Volumes of 
seawater (generally from 1-2 L) are filtered through pre- 
washed, pre- ashed, and pre- weighed 0.7 pm nominal pore 
size GF filters. 

After seawater filtration, the filter (specifically the fil- 
tration area and border) is washed with distilled water and 


stored at — 18°C for subsequent laboratory analysis. Be- 
fore final weighing, the filters are dried at 75°C for 1 h, and 
then temporarily stored in a desiccator. The concentration 
of TSM (in milligrams per liter) is calculated using: 


Ctsm = 


W t - W f 
V w 


w b 


( 20 ) 


where W t is the weight in milligrams of the sample fil- 
ter after filtration; Wf is the weight in milligrams of the 
filter before filtration; and w b is a correction term (in mil- 
ligrams) introduced to account for changes in the weight of 
the filter sample due to changes induced by environmental 
conditions and handling effects between the two weighing 
steps. 

The term w b is obtained from blank filters (i.e., GF 
filters completely conditioned for TSM analysis, not used 
for seawater filtration, but subjected to the same processes 
as the sample filters: transportation to the measurement 
site and back, storage in the freezer, and subsequent filter 
drying). It is computed from the difference between the 
average final weight of the blank filters and their original 
average weight. 

The use of GF filters with a 0.7 pm nominal pore size 
for TSM analysis, can produce underestimates of CtsMi 
because of the loss of particles with diameters less than 
0.7 pm. The filter rinsing for salt removal, plus the fil- 
ter conditioning before and after filtration before weigh- 
ing, however, can induce uncertainties much larger than 
the mass percentage of particles with diameters less than 
0.7 pm. 

An analysis of measurement reproducibility performed 
with 93 pairs of duplicate samples shows average APDs 
equal to 13.9 ± 13.4% for Ctsm— 0.86 ± 0.40 mgL -1 . The 
largest differences (i.e., greater than 30%) between dupli- 
cate samples, have been observed with Ctsm lower than 
0.5mgL _1 and are attributed to the intrinsic uncertainty 
in sample preparation (i.e., water sample nonhomogeneity 
and filter rinsing). 
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Table 9. A summary of the results on the analysis of PSD measurements at different diameter ranges, in 
micrometers, of particles. The average APDs resulting from duplicate analysis of the same sample is d\. The 
values in parentheses indicate the average number of particles N p per cubic centimeter in the given diameter 
range, for the analyzed samples. The standard deviations are indicated by the “d:” sign. 


Diameter Range 

2-3 

3-5 

5-7 

7-10 

di [%] 
n = 230 

17.8 ±18.8 
(20,332 ±12,678) 

9.6 ± 11.1 
(11,346 ±7,588) 

15.8 ± 14.4 
(3,280 ± 2,692) 

25.8 ±22.9 
(1,559 ±1,580) 


Table 10. Chlorophyll a concentration C a) total suspended matter Ctsm> an d colored dissolved organic matter 
absorption coefficient a ys (400) at 1, 8, and 14 m depth from seawater sample analyses for the stations presented 
as case studies in Figs. 7-9. 


Station 

Date 

Depth [m] 

C a [^gL *] 

Ctsm *] 

a ys 

in" 1 ] 

C390401 

26 February 1998 

1 

4.89 

2.33 

0.162 

C390401 

26 February 1998 

8 

1.75 

0.53 

0.106 

C390401 

26 February 1998 

14 

1.51 

0.80 

0.130 

C400401 

10 March 1998 

1 

1.29 

0.53 

0.089 

C400401 

10 March 1998 

8 

1.28 

0.73 

0.106 

C400401 

10 March 1998 

14 

1.38 

1.27 

0.097 

C490201 

9 September 1998 

1 

0.53 

0.27 

0.075 

C490201 

9 September 1998 

8 

0.51 

0.27 

0.013 

C490201 

9 September 1998 

14 

0.91 

0.80 

0.049 


8.10 Particle Size Distributions 

The rationale for making PSD measurements is a desire 
to support theoretical estimates of the seawater particle 
scattering phase function through Mie calculations. The 
PSD (i.e., the number and size of particles suspended in 
seawater) is measured using the Coulter Counter technique 
and making use of a MULTISIZER Ilf system. The Coul- 
ter Counter technique is based on the principle that par- 
ticles, suspended in an electrically conductive liquid and 
forced to flow through a small aperture (orifice) having an 
immersed electrode on each side, change the resistance be- 
tween the electrodes. The latter produces a voltage pulse 
proportional to the particle size, which can be detected 
and measured. 

Analysis of seawater samples was carried out using a 
100 fim aperture, which ensures characterization of parti- 
cles in the diameter range of 1.8-64 ^m distributed over 
256 logarithmically spaced channels (i.e., the diameter as- 
sociated with each channel halves every 50 channels down). 
For each sample, three sequential measurements are car- 
ried out with 2 mL of seawater each time and the resulting 
distribution is averaged. The background of particle size 
distribution is estimated making use of seawater filtered at 
0.22 /im. 

Absolute calibration of the MULTISIZER II system is 
performed annually making use of 5 pm monosized smooth 
particles (i.e., latex formulations of polystyrene divinyl 
benzene) of known density. Measurement repeatability was 
evaluated through duplicate analysis of samples from the 


f Manufactured by Beckman Coulter, Inc., Hialeah, Florida. 


same water volume. The PSD results are given in Table 9 
for different ranges of particle diameters (i.e., 2-3, 3-5, 5- 
7, and 7-10 /mi), and show average APDs ranging from 
9.6-25.8% in the 3-5 ^m and 7-10 /im ranges, respectively. 
The larger average differences observed in the 2-3 /im in- 
terval, when compared to the 3-5 pm interval, are due to 
an increase in the measurement noise at the lower limit of 
the measurement interval, i.e., 1.8 /im. 

9. SAMPLE DATA 

A comprehensive and detailed analysis of the data col- 
lected in the framework of the CoASTS Project for a three- 
year period starting in the latter part of 1995 is given in 
Berthon et al. (2002). For completeness of the overall work, 
an overview of optical and hydrographic profiles typical of 
the AAOT site in the northern Adriatic Sea is given here 
and discussed below. 

In Figs. 7-9, L u (z,\) and c t - w (z, A) at representa- 
tive wavelengths plus S w (z) are displayed for three dif- 
ferent stations: C390401 on 26 February 1998, C400401 
on 10 March 1998, and C490201 on 9 September 1998. 
Relevant variables identifying the seawater biogeochemi- 
cal and optical characteristics of the stations are given in 
Table 10. 

Figure 7 displays the profiles for station C390401 and 
shows some seawater vertical stratification between the 
surface and 3 m depth. Coefficients of c t - w (z,\) show a 
sharp variation of about a factor of four at 412 nm, in a few 
meters range around 3 m depth. In agreement with these 
changes, L u (z, A) profile data (in logarithmic scale) show 
a change in slope at the same depth. The seawater salinity 
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exhibits a layer with approximately 33PSU between the 
surface and 2 m depth, followed by an abrupt increase up 
to 37PSU from 2-6 m. This profile can be explained by 
the existence of a layer of fresh water, rich in particles and 
CDOM originating from the major northern rivers that 
drain into the northern Adriatic Sea. The relatively high 
values of O a , Ctsm> and to a lesser extent of a ys in the 
surface layer (i.e., 1 m depth), given in Table 10 for station 
C390401, clearly support this hypothesis. In addition, the 
high concentration of C a (i.e., 4.89 /igL -1 ) in the same 
surface layer, suggests the presence of nutrients is proba- 
bly due to coastal transport. 

Figure 8 presents the profiles for station C400401, for 
which the water column is vertically unstratified. Both the 
c t - w {z, A) and S w (z) values show very slight changes as a 
function of depth. In full agreement with this, In [L u (z, A)] 
(excluding the data at 665 nm) exhibits a linear decay as 
a function of depth with no significant change in slope. 
With the exception of Ctsm, which shows some slight in- 
crease with depth, the C a and a y8 data for station C400401 
(Table 10) do not exhibit any relevant change with depth, 
which is in complete agreement with the profiles presented 
in Fig. 8. 

Figure 9 displays the profiles for station C490201 and 
shows a stratification near the bottom. The c t - w (z,\) 
data displays almost constant values from the surface up to 
12 m depth, and then exhibits a sharp increase by a factor 
of about 3 from 12-15 m depth. This feature is confirmed 
by a change in the slope of L u (z, A) at 12 m depth. An 
explanation for this could be the existence of a bottom 
current or vertical mixing, because the water is almost 
homogeneous in S w (z ), thereby resuspending the particles. 
The former observations axe supported by the Ctsm data 
at 14 m depth, which exhibit values three times higher than 
those at 1 and 8 m. 

The sample data presented here, although far from 
describing all the possible conditions encountered at the 
AAOT site, give a general overview of the variability that 
can characterize a coastal measurement site in the north- 
ern Adriatic Sea. 

10. CONCLUSIONS 

This report presented the measurement program used 
within the CoASTS Project (at the AAOT site) and exten- 
sively discussed the applied methods with a particular em- 
phasis on those more relevant for bio-optical modeling and 
the calibration and validation of ocean color sensors. The 
methods discussed for the analysis of radiometric measure- 
ments highlighted the importance of correcting in-water 
radiance and irradiance data from superstructure (tower- 
shading) effects, instrument self-shading, and bottom ef- 
fects. These perturbations for the specific measurement 
conditions analyzed in this study, may induce an overall 
spectral uncertainty in L u (0"A) ranging from 4-18% at 
412 nm, —2 to 4-10% at 490 nm, —9 to +10% at 555 nm, 
and 6-16% at 665 nm. 


The intercomparison of seawater absorption coefficients 
from AC-9 in situ measurements and from laboratory spec- 
trometric analysis of discrete seawater samples, showed the 
best agreement when the scattering correction— for per- 
turbation effects induced in AC-9 data by the finite ac- 
ceptance angle of the optics and the non-ideal reflective 
surface of the absorption chamber — is performed removing 
a variable percentage of the scattering coefficient. In this 
case, the scatterplots of the absorption coefficients from 
spectrometric analysis of seawater samples and AC-9 mea- 
surements showed average absolute percentage differences 
of _2, -19, -18, and -18% at 412, 488, 555, and 676 nm, 
respectively. 

The intercomparison of the spectrometric transmission 
method with the T-R method used in CoASTS for the 
analysis of seawater particulate absorption a p (A), showed 
differences ranging from +11% up to —15% in the 412- 
665 nm spectral interval. The reproducibility of CoASTS 
particulate and yellow substance absorption coefficients, 
determined with an analysis of replicate samples, showed 
average APDs of 9.5 + 5.5% for a p (412) = 0.093+0.051 m" 1 
and 12.1 + 6.3% for a y *(412) = 0.175 + 0.038 m“\ respec- 
tively. 

A measurement reproducibility study, using replicate 
samples, was also presented in determining the concen- 
tration of the most represented pigments in the northern 
Adriatic Sea obtained with HPLC analysis. The results 
for C a showed average APDs of 7.1 ± 6.3% for C a =0.86± 
0.19 ^gL -1 . A similar analysis was carried out to assess 
the reproducibility of Ctsm measurements. The latter re- 
sults showed average APDs of 13.9+13.4% for Ctsm— 
0.86 ± 0.40 mgL“ x , exhibiting the largest differences (i.e., 
up to 30%) for concentrations lower than 0.5 mgL -1 . 

The study on measurement reproducibility highlighted 
a strong variability for all quantities obtained from the 
analysis of replicate samples produced from the same vol- 
ume of seawater, and suggested difficulties in properly ho- 
mogenizing the water volume during the preparation of the 
samples. 

The uncertainty, repeatability, and reproducibility anal- 
yses presented in this study, were a first attempt to esti- 
mate the accuracies in the different quantities measured 
within the CoASTS Project. Although some of the results 
may suggest an unsatisfactory accuracy for specific appli- 
cations, it must be remembered that the coastal nature 
of the study area significantly increases the measurement 
difficulties, and in many cases, may reduce the achievable 
accuracies. 

The sample profile data, gave an overview of the typi- 
cal seawater bio-optical characteristics at the measurement 
site. Their variability, assumed to be quite representative 
of the whole northern Adriatic Sea, confirms the possibility 
of exploring very different water types at a single measure- 
ment point and, thus, sustains the suitability of the AAOT 
measurements to support bio-optical modeling, as well as 
calibration and validation activities for northern Adriatic 
Sea coastal waters. 
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Glossary 

AAOT Acqua Alta Oceanographic Tower 
A/D Analog- to- Digital 

ADEOS Advanced Earth Observing Satellite 
AERONET Aerosol Robotic Network 

APD Absolute Percent Difference 

ARGOS Not an acronym, but the name given to the 
data collection and location system on the 
NOAA operational satellites. 

BOUSSOLE Boude pour V acquisition de Series Optiques a 
Long Terme (buoy for the acquisition of a long- 
term optical series). 

CDOM Colored Dissolved Organic Matter 
CNR Consiglio Nazionale delle Ricerche (Italy) 

CTD Conductivity, Temperature, and Depth 
CoASTS Coastal Atmosphere and Sea Time Series 
CZCS Coastal Zone Color Scanner 

DCP Data Collection Platform 
DHI DHI Water and Environment Institute (Den- 
mark) 

FEL Not an acronym, but a lamp designator. 

GF Glass Fiber Filter 
GSFC Goddard Space Flight Center 

HPLC High Performance Liquid Chromatography 
IS Internal Standard 

ISDGM Istituto per lo Studio della Dinamica delle 
Grandi Masse (Italy) 

JGOFS Joint Global Ocean Flux Study 
JRC Joint Research Centre (Italy) 

LOV Laboratoire d’Octanographie de Villefranche 
(France) 

LPCM Laboratoire de Physique et Chimie Marines 
(France) 

LUT Look-Up Table 

MERIS Medium Resolution Imaging Spectrometer 
METEOSAT Meteorological Satellite 

MOBY Marine Optical Buoy 


MODIS Moderate Resolution Imaging Spectroradiome- 
ter 

MOS Modular Opto-electronic Sensor 

MVDS Multichannel Visible Detector System 

NASA National Aeronautics and Space Administra- 
tion 

NAS DA National Space Development Agency (Japan) 
NIST National Institute of Standards and Technol- 
ogy 

NOAA National Oceanographic and Atmospheric Ad- 
ministration 

OCI Ocean Color Irradiance 
OCR Ocean Color Radiance 

OCTS Ocean Color and Temperature Scanner 

PlyMBODy Plymouth Marine Bio-Optical Data Buoy 
PML Plymouth Marine Laboratory (England) 

POLDER Polarization and Directionality of the Earth’s 
Reflectance 

PSD Particle Size Distribution 

SeaWiFS Sea-viewing Wide Field-of-view Sensor 
SNR Signal-to-Noise Ratio 

T Transmission method for spectrophotometric 
analysis. 

T-R Transmission- Reflection (method for spectro- 
photometric analysis) 

TSM Total Suspended Matter 

VKI VKI Institute for Water Environment (Den- 
mark) 

WETLabs Western Environmental Technology Laborato- 
ries (Inc.) 

WiSPER Wire-Stabilized Profiling Environmental Radi- 
ometer 

YBOM Yamato Bank Optical Mooring (Japan) 

YES Yankee Environmental Systems (Inc.) 

Symbols 

Indication of a measured quantity with any giv- 
en variable. 

a(z, A) Total seawater absorption coefficient at depth 
z. 

Qdp( A) Nonpigmented particulate matter absorption 
coefficient. 

Qdp{z, A) Nonpigmented particulate matter absorption 
coefficient at depth z. 

a p (A) Particulate matter absorption coefficient. 

a p (z y A) Particulate matter absorption coefficient at 
depth 2 . 

a ph(A) Pigmented particulate matter absorption coef- 
ficient. 

d p h(z y A) Pigmented particulate matter absorption coef- 
ficient at depth z. 

at- w (z , A) a^J U) (z, A) corrected for AC-9 scattering effects. 
cit-w(z, A) Total seawater minus pure seawater absorption 

^ coefficient from AC-9. 

&t-w(z } A) afT w (z y A) corrected for salinity and tempera- 
ture effects. 

a pjJ (A) Colored dissolved organic matter absorption co- 
efficient. 

a ys ( 2 , A) Colored dissolved organic matter absorption co- 
efficient at depth 2 . 
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a™ (A) Pure seawater beam absorption coefficient. 
av Average. 

A Intercept. 

A p Pigment peak area. 

A r Reference internal standard peak area. 

A s Sample internal standard peak area. 

A s1iS (A) Equivalent particle suspension absorbance. 

Ays (A) Dissolved matter absorbance. 

6(A) Total seawater scattering coefficient. 
b t -w(z, A) Total minus pure seawater scattering coefficient. 
B Slope. 

B e Buffer-extract dilution factor. 
c(z, A) Total seawater beam attenuation coefficient. 
c t -w{z, A) Total seawater minus pure seawater beam at- 
tenuation coefficient. 

ct- w {z, A) Total seawater minus pure seawater beam at- 
tenuation coefficient from AC-9. 
cf-w(z, A) c t -xv(z, A) corrected for salinity and tempera- 
ture effects. 

Cu,(A) Pure seawater beam attenuation coefficient. 

C a Chlorophyll a concentration. 

Cg d ( A) Absolute calibration coefficient for E d (Zy A). 

Cg (A) Absolute calibration coefficient for E u (z, A). 

Ci (A) Absolute calibration coefficient for L u (z, A). 

C p Pigment concentration. 

C^(A) Absolute calibration coefficient for 3( 2 4)- 
Ctsm The concentration of total suspended matter. 
CC Cloud cover. 

d\ Average differences from duplicate analysis of 
the same sample. 

d 2 Average differences from analysis of duplicate 
samples. 

da Average differences from duplicate samples an- 
alyzed by different laboratories. 

D^(X) Average dark value for V^(z y A). 

D y Sequential day of the year. 

E{ A) Direct sun irradiance. 

Eo(X) Extra-atmospheric sun irradiance. 

Ed{z, A) Downward irradiance at depth z. 

E d (z, A) E d {z y A) affected by measurement perturba- 
tions. 

E d {0~,X) Subsurface downward irradiance. 

Ed{ 0 + ,A) Above- water downward irradiance. 

Ed( 0”,A) E d (0~ y \) affected by measurement perturba- 
tions. 

E t ( 0 + ,A) Above-water downward diffuse irradiance. 

E s ( A) Direct normal sun irradiance. 

E u {z y A) Upward irradiance at depth z. 

E u (z, A) E u (z, A) affected by measurements perturba- 
tion. 

E u ( 0", A) Subsurface upward irradiance. 

E u (Q-,\) E u ( O’, A) affected by measurement perturba- 
tions. 

Eu C { 0~>A) Bottom corrected upwelling irradiance. 
f p Pigment specific calibration factor. 

F a Filter clearance area. 

Lj (A) Immersion coefficient. 

3(z,A) Absolute calibration of the considered radio- 
metric quantities. 

3(0", A) Subsurface radiometric quantity. 


3(0", A) Subsurface radiometric quantity affected by 
measurement perturbations. 

3(0, A) Radiometric quantity affected by measurement 
perturbations. 

A:(zb, A) Transmittance for the downward plus upward 
normal optical paths beween the surface and 
the bottom at depth z. 

KE d { A) Diffuse attenuation coefficient from E d (ZyX). 

K A) Diffuse attenuation coefficient from Ed{z } A). 
K.E d { A) Average profile diffuse attenuation coefficient 
from Ed{z, A). 

K % (A) Diffuse attenuation coefficient from E u (z, A). 

K^{ A) Diffuse attenuation coefficient from 3(^4). 

K l (A) Diffuse attenuation coefficient from L u (z,A). 
L c Cuvette pathlength. 

Li(0, 4 > , A) Diffuse sky radiance. 

L u {z,\) Up welled radiance at depth z. 

L u (z, A) L u (z, A) affected by measurement perturba- 
tions. 

Lu(0‘,A) Subsurface upwelled radiance. 

Tti(0",A) L u ( 0",A) affected by measurement perturba- 
tions. 

L u ( 0 + , A) Above- water upwelled radiance. 

Lu°( 0”, A) Bottom corrected upwelled radiance. 

Lwn( A) Normalized water-leaving radiance. 
m The relative air mass. 

M Sea state, 
n Number of points. 

N Number of optical stations. 

N p Number of particles. 

P a Atmospheric pressure. 

Q n {z, A) E u (z, A)/L u (z, A),i.e., the Q-factor at nadir. 
Q n (0',A) E U (Q~ y A)/ L u (0~ , A). 

$40", A) E u (0~ y A)/ Lu(0“, A). 

r 2 (A) Determination coefficient. 

r*(A) Ratio between diffuse and direct above- water 
downward irradiance. 

R(zy\) Bottom correction of irradiance reflectance. 

R( 0",A) E u {Q'yX)/E d (0~yX). 

Ri Instrument case radius. 

Rs Sensor entrance optics radius. 

RH Relative humidity. 

S w (z) Seawater salinity at depth z. 

SD Secchi disk depth. 

T a Air temperature. 

Ti Tide level. 

T w (z) Water temperature at depth z. 

V e Volume of pigment extract. 

Vi HPLC injected volume. 

Vj(z, A) The 3(z, A) voltage (in digital counts). 

V w Volume of filtered water. 

Wb Correction term. 

W d Wind direction. 

Wf Filter weight before filtration. 

W a Wind speed. 

Wt Filter weight after seawater filtration, 
z Depth. 

zb Bottom depth. 

zo Upper limit of the extrapolation depth. 
z\ Lower limit of the extrapolation depth. 
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5 Fixed percentage. 

773 (A) Perturbation correction factor for 3(CT,A). 

Bottom effects correction factor for 3(0~, A). 
r/^ 5 (A) Self-shading correction factor for 0(0”, A). 
t?J 5 (A) Tower-shading correction factor for 3(0”, A). 

0 Zenith angle. 

9q Sun zenith angle. 

A Wavelength. 

Ao Reference wavelength. 

Pb{ A) Bottom reflectance. 

a Standard deviation. 

r Atmospheric optical thickness. 
r a (A) Aerosol optical thickness. 
t 0 (A) Ozone optical thickness. 
t;?(A) Rayleigh (molecular) optical thickness. 

4> Azimuth angle. 

0o Solar azimuth angle. 

t^o(A) Single scattering albedo. 
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